Generalized Quantum Baker Maps as perturbations of a simple kernel 
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We present a broad family of quantum baker maps that generahze the proposal of Schack and Caves to any 
even Hilbert space with arbitrary boundary conditions. We identify a structure, common to all maps consisting 
of a simple kernel perturbed by diffraction effects. This "essential" baker's map has a different semiclassical 
limit and can be diagonalized analytically for Hilbert spaces spanned by qubits. In all cases this kernel pro- 
vides an accurate approximation to the spectral properties - eigenvalues and eigenfunctions - of all the different 
quantizations. 



PACS numbers: 

I. INTRODUCTION 

The accurate description of pure eigenstates of chaotic sys- 
tems is still beyond us. The initial Berry- Voros ansatz 1 1] as- 
suming a microcanonical "uniformity" on the energy shell in 
accordance with Schmirelman theorem leaves space for 
weaker structures that lead to scaiTing by unstable periodic or- 
bits 01 ■ The choice of basis is of course crucial in achieving a 
simple description. The standard choice for numerical calcu- 
lations is to diagonalize the hamiltonian - or the map - in an in- 
tegrable basis: coordinate, momentum, oscillator, plane wave, 
etc. In these bases the chaotic eigenfunctions have very large 
participation ratios, of the order of the Hilbert space dimen- 
sion, assimilating them to random eigenfunctions and exclud- 
ing any perturbative description. Another approach is to use 
a basis of "quasi-modes" constructed on short periodic orbits 
, which has been used extensively to describe the eigen- 
states of the Bunimovich stadium. Still another approach has 
been the study of the distribution of zeroes in analytic repre- 
sentations, p*! 

In this paper we show that the eigenfunctions and eigenval- 
ues of the quantized baker's map - and a great variety of its 
generalizations - can be approximated by a relatively simple 
basis which diagonalizes a kind of "essential" baker's map. In 
this basis the eigenfunctions of all quantizations of the baker's 
map have very small participation ratios indicating that the ba- 
sis vectors are excellent approximations to the eigenfunctions. 
Moreover the eigenvalues are well reproduced in first order 
perturbation theory. This basis can be analytically constructed 
for special Hilbert space dimensions, notably the qubit case 
when D ^ 2^ |6]. Our approach extends the pioneering 
work by Lakshminarayan |7, 8], who was the first to realize 
that some eigenfunctions had a simple structure when looked 
upon differently in the Walsh-Hadamard basis. 

The organization of the paper is as follows: in Section II 
we review the quantization of the baker's map and propose 
a wide generalization of the Schack and Caves [9] construc- 



tion for qubits by allowing arbitrary Hilbert space dimensions 
and Floquet angles. Utilizing a mixture of analytical and cir- 
cuit techniques we give a semiclassical interpretation of this 
construction in terms of semiquantum maps that quantize the 
baker's map iterates 1 10]. We show that all these maps share 
a common structure as the product of two unitary kernels: a 
fixed one *B common to all families and an almost diagonal 
one containing diffraction effects. In Sec. Ill we analyze spec- 
tral properties of this family of maps and show that S pro- 
vides an accurate basis for the description of eigenfunctions 
and eigenvalues. When the dimension is D = 2^ this ba- 
sis can be analytically constructed and is labeled by primitive 
binary strings. For maps with antisymmetric boundary condi- 
tions this basis improve the Hadamard representation discov- 
ered by Lakshminarayan |^7J. 



II. BAKER'S MAP QUANTIZATIONS 
A. Quantum baker maps 

The classical baker map is an example of an intuitive canon- 
ical transformation that can be expressed in terms of symbolic 
dynamics using the binary Bernoulli shift. The map is defined 
in the unit square phase space {q,p e [0, 1)) by 

= 2q, - [2q,] (1) 
K+l = {p, + [2q,])/2 (2) 

where [ ] denotes the integer part, and i is the discrete time. 
This map is area-preserving, and geometrically stretches the 
square by a factor of two in the q direction, squeezes by a 
factor of a half in the p direction, and then stacks the right 
half onto the left. 

The map has a simple symbolic dynamics involving the 
binary expansions of the coordinates, q ~ O.eoei ... = 
EZo^k2-'-' andp = 0.e_ie_2... = EZi^-k^'' 
{e-i e 0, 1). One point of the phase space is represented by 
a bi-infinite symbolic string 
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(p,q) ^ . .. e-2e-i • eoeie2e3 



(3) 



2 



and the baker map action upon symbols is 

{p,q) — > {p',q') = 



e-2e-ieo • eie2e3 ■ 



(4) 



The classical map has two symmetries: time-reversal (T) 
reversing the direction of the flow and exchanging p ^ q; 
and parity (R) exchanging q ^ 1 — q,p ^ 1— and bitwise 
logical NOT (0 ^ 1) upon symbols. 

A quantization of the baker map is well known and pro- 
ceeds by first setting up a Hilbert space with appropriate 
boundary conditions. The phase space is made compact by 
imposing perodic boundary conditions, thus turning it into a 
2-torus. The corresponding quantum structure is character- 
ized by quasi-periodic boundary conditions 

^{q+l) = e'^^yj{q), i>{p + 1) = e~'^^^^{p) (5) 

where 2TTri and 2ttk are arbitrary Floquet angles and are 
Fourier transformed pairs. Solutions to (|5} only exist if hD — 
1 with D integer, and they span a Hilbert space H]^'^ having 
finite dimension D. The position and momentum eingevectors 
are \qj) and \pj) with eigenvalues qj = [j + rj)/D andpj. 
{k + k)/D respectively (j, /c — 0, . . . ,D — 1). The vectors 
of each basis are orthonormal, {qj\qk) — {Pj\Pk) = Sjk and 
the two bases are related via the finite Fourier transform with 
arbitrary Floquet angles. 



(6) 



The quantization of the baker map on an even-dimensional 
Hilbert space can be achieved 1 11 ] converting the most signif- 
icant bit of position in the most significant bit of momentum. 
The matrix of the map in mixed representation has the form 
of two blocks with the finite Fourier transform of size D/2m 
each one. In position representation we have 



^pos — 




(7) 



This matrix product has a simple circuit representation in 
terms of the Fourier transform. This is shown in (Fig. [Q, 
where the lines represent subspaces ordered with the most sig- 
nificant one on the bottom, the box is a unitary operator acting 
in the respective space, and the temporal flux is from left to 
right (opposite to the matrix representation) 1 12]. 

The classical symmetries are preserved for some values of 
the Floquet angles. T is preserved when r/ = k, and R when 
T] + K — 1. The original quantization with rj = k — was 
made by Balazs and Voros fTHl . In particular, the choice rj — 
K = 1/2 (antisymmetric boundary conditions) yields a map 
with both symmetries 1 13]. 

We will use throughout techniques borrowed from quantum 
information theory ] 12] allowing simple graphical representa- 
tions of unitary operators in tensor product Hilbert spaces. For 
the baker's map the well known decomposition of the Fourier 
matrix into qubit operations ] 14] when D — 2^ allows a sim- 
ple and efficient representation in terms of circuits. In Ap- 
pendix A we show that this decomposition can be extended to 
arbitrary Floquet angles and other factorizations of D. 
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Figure 1: Circuit representation of the Balazs-Voros-Saraceno quan- 
tum baker map. The thick line represents a subsystem of dimension 
D/2. 



B. Family of baker maps on qubits 

The original quantization scheme only required D to be 
even. In the special case when D ~ 2^ , Schack and Caves 

fjposed an entire class of quantum baker maps on qubits 
[l51 . They connected the binary representation of the clas- 
sical baker map to the qubit structure using the partial Fourier 
transform, defined in the general case as 



G!?'" EE U 



n = 1, . 



(8) 



where 12" is the unit operator on the first n qubits, and F^lj^_„ 
is the Fourier transform on the remaining qubits. This trans- 
formation is used to define orthonormal basis states that are 
localized on the unit-square phase space [ISJ. 

This operator has a straightforward representation in terms 
of quantum circuits as shown in Fig. 2. 



N-n qubits 



n qubits 



F^N-n 



Figure 2: Circuit representation of the partial Fourier transform act- 
ing on the N-n less significant qubits. 

A quantum baker map with N qubits is defined for each 
value of n = 1 , . . . , iV by the unitary operator |^ 



B 



N,; 



Gn-1 O Sn O G„ 



(9) 



where the dependence in rj and k is implicit. The 
shift operator Sn acts cyclically only on the first 
n qubits, i.e., S'„|xi)|a;2) . . . |a;„)|a::„+i) . . . |a;Ar) = 
12^2) • ■ • |a::„)|a;i)|a::„+i) . . . \xn)- 

Since is the unit operator, B^ i = Bbvs is the Balazs- 
Voros-Saraceno (BVS) map. On the other hand, it is easy to 
show that Bm.n is a map constructed only with swaps and 
one qubit Fourier transform. Since Sn commutes with Gn^, 
Bi\f ,n, it can be written as 



-Bat „ — 



B 



N-n+1,1 



O Sn- 



(10) 



3 



Thus the action of ijjv.n is equivalent to a cyclic shift of the 
n most significant qubits followed by application of Bbvs to 
the N — n + 1 least significant qubits. This equivalence is 
shown in circuit representation in Figs. |3]and0] All these 

















Figure 3: Circuit representation of BN,n in terms of the BVS baker 
map 

maps have efficient implementations in terms of single qubit 
operations, swaps and controlled phases, even for arbitrary 
Floquet angles, (see Appendix lAt. 
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C. Quantization of iterated map 

We can give a different -semiclassical- interpretation to the 
Schack and Caves construction that, besides clearly showing 
that all maps in the family are equivalent in the semiclassical 
limit (h ^ ^ 0), allows a generalization to Hilbert space 
dimensions of a much wider class. 

We start by recalling that the scheme utilized by Balazs and 
Voros inil to quantize the original baker map can also be uti- 
lized to quantize its iterates. We denote B^"^^ the unitary map 
resulting from the direct quantization of the T iteration of the 
classical map. Since the classical map divides the q or p seg- 
ments into 2'^ equal intervals, we should initially assume that 
divides D so that each strip is quantized by states. 
Notice that for this construction to be semiclassically plausi- 
ble there should be many quantum states in each strip thus 
requiring 00. As an example, two iterations of the 

classical baker map require four strips as in Fig. [S] 

V, V, V, V. 
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Figure 4: Circuit representation of the four members of the quantum 
baker map for four qubits {q = k = G). H is the Hadamard operator 



We will dedicate ourselves to B^^n, the last map of this 
family, in section D. This map has appeared in different con- 
texts. In the literature on quantum walks fl^ it is called the 
"many coins" map and is used to study quantum walks with 
many independent coin throws. In the context of quantum 
graphs it represents a De Bruijn graph with complex phases on 
its edges flT]. Nonnenmacher has observed that if the Fourier 
transform is replaced by the Walsh-Hadamard transform in 
the Balasz-Voros scheme then B^.n is the Walsh quantiza- 
tion of the baker's map and constitutes a "toy model" for the 
baker's map 1 18]. Taking this observation into account we can 
also think of the intermediate members of the family as mixed 
Fourier - Walsh-Hadamard quantizations. 

If one is willing to admit a BVS map in dimension D ~ 1, 
which would amount to a trivial inessential phase, then Sn - 
the simple cyclic shift on iV-qubits - could also be considered 
as an extreme member of this family. It was in fact a quanti- 
zation of the baker's map proposed by Penrose fl^ . but never 
pursued. 



Figure 5: Representation in phase space of two iterations of the 
classical baker map 

The quantization of this second iterate proceeds exactly as 
for the first. The map B^^'' can be quantized in the mixed 
representation as 



B. 



r^,K (2) 
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(11) 



These quantizations for various values of T were studied 
in 1 10|, and provide a semiquantum dynamics, somewhat in- 
termediate between semiclassical and quantum propagation. 
The matrix representation of B^^^ can be expressed easily us- 
ing finite strings of bits in symbolic dynamics. 



{m\Bi 



(12) 

where u{qo) = (eo . . . er-i) = and u''{qo) = (et-i • ■ • eo)- 

The tensor product sctructure of i?^^) allows a very simple 
circuit representation (Fig. |6|l, using the swap operator 82- 
The inversion of the order of the most significant qubits cor- 
responds to the transfer of a finite string of T bits of symbolic 
dynamics from momentum to coordinate and is accomplished 
by ^(2^) swaps for even (odd) T. 

It is important to observe that the direct quantization of the 
T-step classical map and the T-th iteration of the one-step 
quantum map do not produce exactly the same matrices. In 
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It should be clear that many more inequivalent construc- 
tions are posible along these lines, i.e. ij(3)ij(i)t5(i)t ^nd 
many other similar combinations. 

This interpretation of the family as a product of semiquan- 
tum maps gives also a natural explanation to the result of 1 2^ 
who found that Bn^n was abnormal in that its classical limit 
was different from the baker map. Seen from a semiclassical 
point of view, this result is quite natural as it would involve 
the quantization of strips of dimension = 1 by means 

of rank one projection operators, thus strongly violating any 
semiclassical justification. However, as we show in Sec. Ill, 
it is precisely this "extreme" quantum map that provides the 
approximate description of all the other quantizations. 



Figure 6: Circuits in the mixed representation for the quantization 
of the first four iterations of the baker's map. 



fact quantization and propagation are operations whose com- 
mutativity is only justified in the semiclassical limit D — > 
oo. The Schack and Caves construction can then be seen 
from a different perspective using the fact that ) and 

B^B^ represent different matrices which are only asymptot- 
ically equal whenever the semiclassical condition is satisfied. 
Thus we can define families of quantizations of the one step 
map as 



B 



J5(n-l)t^(n) 



l,...iV 



(13) 



where the only requisite is that 2^ divides D, and the BVS 
map is recovered for I3d,i- It is clear from the semiclassical 
arguments in jToll that all these families are semiclassically 
equivalent in the sense that they all represent valid, and differ- 
ent, quantizations of the one-step baker map. However, this 
equivalence requires that n remain fixed while ^ ^ oo. If, 
on the other hand we allow 2" ^ Z) we can expect strong 
quantum deviations from this equivalence. 

If we work with a Hilbert space spanned by N qubits and 
for antiperiodic boundary conditions, it is easy to see that the 
quantization obtained as a product of two semiquantum maps 
is equivalent to the family Bat „ of Eq|9] obtained by Schack 
and Caves. In Fig. for example, we use circuits to show 
thatE(3)tB(4) ^ Bn4. 
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Figure 7: Equivalence of the quantum baker map with the semiclas- 
sical definition and with the Schack and Caves maps. 

In what follows we will use arbitrary Hilbert space dimen- 
sions. Therefore we will label the different maps not by the 
number of qubits but by this dimension. Thus the Schack and 
Caves family will be labeled as _B2« 2'>-i. 



D. Quantum baker's map generalization 

Motivated by the Schack and Caves families for a space 
spanned by qubits (Eq. |9}, we will generalize the quan- 
tum baker's map (QBM) to Hilbert spaces with arbitrary Flo- 
quet angles and arbitrary even dimension. Assume that D = 
2DaDp in which case we interpret the map as operating in a 
product Hilbert space Tin — 0.2 ^Ti-n^ ^Hd^- We relabel 
the coordinate states as 



\j) |ej/3ia) = |e) ® li/3) \ja) 



(14) 



where j = eDaDp + jpDa + ja and e = 0, 1; < j/j < Di3\ 
1^ ja < Da thus making 7^2 the "most significant" sub- 
space and Ha the least significant one. The QBM families are 
defined as 



B 



D,Dp 



= 1 



{idp®b) 



(15) 



where B is the BVS baker's map in a 2Z3q -dimensional 
Hilbert space with implicit dependence on the Floquet angles. 
The shift 52, d,, between a qubit and a j3 subspace is a partic- 
ular case of a shift defined as 



S 



D-l 

D1.D2 = X! 

J=0 



D^j-[^]{D^\)) 01 (16) 



where [j] is the integer part of j. This shift is a permutation 
of the states j £ Hd which exchanges the significance of the 
Di, D2 subspaces. 

Bo.Dp has a very simple circuit representation (Fig. |8|i 
which obviously generalizes the qubit circuit of Fig. 
Clearly different maps with a constant value of DaDp con- 
stitute a family of quantizations that generalize the previous 
qubit construction, which is recovered by choosing antiperi- 
odic boundary conditions and a Hilbert space spanned by 
qubits. 

We now arrive at the main point of our construction: the 
baker's map in Fig. (Sjcan be split, using the factorization 
properties of the discrete Fourier transform (cf. Fig. 1171 in 
Appendix A) as the product of two unitary kernels (Fig|9jl. 
The first one which can be considered as an "essential" baker's 
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Figure 8: Circuit representation for the quantum baker map families. 
Each line represents a subspace with the most significant one at the 
bottom. 
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Figure 9: Circuit representation for the QBM family defined in Eq. 
1151 Each line represents a subspace and the squares respresent the 
phase interaction (see App. A). 



map, is an obviuos generalization of the "many coins" map of 
Fig. |4] and is common to all families. We define it as 



*B = e-'^^'^^ So (f^^Id^d, 



-i2-KTjK, { ^ 



D^Db 



where 



is the shift operator used in 0,0]. 

The position representation of is a 

2DaDp X 2DaD/j complex matrix. 
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dj 



(17a) 
(17b) 

(18) 
square 

(19) 



where 



a b 
c d 



V2 



e-"''- I 1 e*""* 

„iTr?7 i7r(j)+K+l~ 



(20) 



The second kernel in Fig. |9]contains all the peculiarities of 
the different family members. It is diagonal in the subspace 
Tip, while the interaction of the qubit line with subspace Tia 
gives raise to diffraction effects. The matrix elements of this 
kernel are 



where 



with 



.v^^ sin(7rl)a7) 
Da sin (7r7) 



(22a) 



7 



Da 



(22b) 
(22c) 



The strong forward diffraction peak gives this kernel an al- 
most diagonal structure with weak off diagonal elements re- 
flecting the block structure of the different sizes of subspaces 
Ha, Ti-p. When Da = 1, the diffraction kernel is the identity, 
and for the qubit case this is the "many coins" map lll6ll2lll . 

The generalizations of the QBM can be rewritten in terms 
of both kernels defined in Eas. ll7all21l andl22las 



B 



D,Df3 



= (Id, ®K)o^ 



(23) 



where -B^i i = B is the Bbvs', Bjy jj /2 = ^ is the essen- 
tial quantum baker's map (EQBM). A similar decomposition 
was obtained in |2^ for the qubit case. 

The modulus of the matrix elements of both kernels for 
Da = 5, Dp = 3 are shown in Fig. 1101 





Figure 10: Modulus of the matrix elements of the diffraction kernel 
(left) and EQBM (right). D = 2D„D^ with Da ^ 5, D/s = S and 
the Floquet angles 7; = k = 0.5 

It should be noticed that, even when antiperiodic boundary 
conditions are chosen, this family of quantizations does not 
respect the time reversal symmetry of the original Bbvs 
map. On the other hand they do commute with the parity R. 

An alternative splitting of the QBM family can be consid- 
ered if the i^2~^ gate in Fig. [Ojis included in the diffraction 
kernel. In this case the resulting "essential" baker map is the 
shift map considered in (tIi. We discuss the relative merits of 
these two splittings in the conclusions. 



III. SPECTRAL PROPERTIES 



(21) 



Based on the fact that the diffraction kernel is almost di- 
agonal we can attempt to describe the spectral properties of 
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the family Sd^d^j in terms of those of the simpler 03. That 
this is not a hopeless endeavour is shown in Fig. n]where we 
compare the positive parity spectra of the family with D = A8 
and various values of D^, Dp. All quantizations share visu- 
ally similar gaps and fluctuations indicating that a common 
structure is present. This fact is even more apparent when the 
spectra are smoothed. The resulting spectral density is almost 
identical for all the maps. 
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Figure 11: Comparison of the positive parity spectra of the fam- 
ily Bn.Dff for D = 48. The values of {Da, Dp) are on the right. 
The top spectrum corresponds to 25 and the bottom to Bvbs- The 
continuous lines are the smoothed spectral densities. 

A more quantitative measure of the similarity of the eigen- 
functions is given by the average participation ratio. The par- 
ticipation ratio, PR — | {fr\'4') l**)^^, is a rough measure 
of the number of basis elements \ifr) needed to construct the 
state with PR e [1, D\. To assess the overall complexity 
of the eigenstates in a given basis we compute the 
average PR 



1 D-l /D-l \ 
a=0 \ i=0 / 



(24) 



This quantity is plotted in Fig. El for Bbvs eigenstates 
(|$Q,)) using several bases and varying the dimension of the 
Hilbert space. 

The (PR) in the position basis grows according to the ran- 
dom matrix theory predictions (just below D/2) revealing the 
chaotic nature of the QBM. The (PR) decreases significantly 
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Figure 12: (Color online) Participation ratio average as a function of 
the dimension for Bbvs with antiperiodic boundary conditions in 
different bases: position basis, shift basis, symmetrized shift basis 
and EQBM basis. 



using the shift eigenstates as a basis ff\, but has strong depen- 
dences on the dimension of the map which can be associated 
to periodic orbits of the map and number theory. The shift 

■■ X^TJq^ e~'^^''^i'j) where are 



eigenstates are \Xj^ 



1 



the trajectory points £)/2l^o)- 

The symmetrized shift, used in L^, is constructed as the 
parity projection of the shift basis. 



1 



(25) 



where i? I j') = jZ?— 1 —j) is the parity acting in position basis. 

This basis has a lower value of the (PR) and less fluctu- 
ations with the dimension D. In the case of EQBM's eigen- 
states as a basis, the (PR) is much lower and grows slowly 
and smoothly with D. The fact that (PR) is extremely small 
for the EQBM basis shows that it should be possible to de- 
scribe the QBM spectral properties (eigenfunctions and eigen- 
values) in terms of the EQBM basis. 

We have also explored the dependence of the growth of 
(PR) in the EQBM basis for different values of the Floquet 
parameters ((77, k)). Values close to the antisymmetric case 
(1/2, 1/2) yield consistently smaller values. 



When D 
in Refs. E 



A. The special case D = 2 . 

2^ the EQBM is the "toy model" considered 
l23il and its spectral properties can be analyzed 
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exactly. Beginning with the shghtly more general operator 
t/ = (m (g) I2JV-1) oS = S o (l2iv-i (E) u) where S is the shift 
operator defined in Eq. ^]and ii is any unitary operator in a 
qubit space. The unitary matrix u can be diagonalized as 





A 



At 



(26) 



where A is a unitary 2x2 matrix. thus, the uni- 
tary matrix U can be written in a simpler way as 
U — A^'^UqA''^^ with Ua in computational matrix repre- 
sentation as 



(ji|C/o|jo> 



exp iUq 



Jo 



2N-1 



xS ji - 2jo + 



30 



2N-1 



(2^ - 1) 1(27) 



This operator takes a very simple form if we label the 
states by a binary string, \j) — > Iv) = |ao...aAr_i) with 



j — ^f^Q^cLi2^ ^ In that representation 



giao6»i+i(l-ao)6»o 1^^ 



■ aN-icioj 
(28) 

This is a permutation which will have cycles whose period T77 
will depend on the binary structure of v. We then have 

U^^\v) = e'^^\v) (29) 

where $77 = NqOq + NiOi, T-jj is the primitive period of the 
string V and A'o and iVi are the number of O's and I's in the 
primitive string. Notice that A^'o + A^i = IV and that T— must 
divide N . The eigenvalues of C/q are then the shifted roots of 
unity 



(30) 



and the corresponding eigenstates are 



Ptl. ^ 



m=0 



where 



$- + 27rfc 



m — 1 

1=0 



(31a) 



(31b) 



For example the primitive orbits for N ~ 4 are 0, 1, 01, 
0011, 0001 and Gill leading to two fixed points, one cycle of 
period 2 and 3 cycles of period 4. 

The eigenstates of the original operator U are then obtained 



(32) 



When u = e '27r»)K^t ('?■«;) ^^.^ eigenvalues and eigenfunc- 
tions are given explicitly as 

i?o,i e*''°'i =±exp(^i|(?7 + K-2?7K)=F«w) (33) 

= — 2=F\/2cos(^|(?7 + k)=f 
e^'^lO) + {±V2e'i^'''+'''>^"^ - 1)|1) (34) 



where w €! [0, 7r/4] is defined as 



V2 



(35) 



These expressions then provide the spectral properties of the 
EQBM in the qubit case D = 2^ for arbitrary Floquet angles. 
For special values of 77, k the matrix -F2 has a definite 

t("'°)^' = land (fV^-'-'-'A' = 1 



short period, i.e. (^Fj y — ^ anu 2 
thus leading to a periodicity 2N or AN for the full EQBM 
map. These short periodicities imply that the spectrum given 
by Eq. |30|will be highly degenerate for large values of N . 
Therefore the eigenfunctions in Eq. |^will not be unique 
and other linear superpositions can also be constructed. They 
have been considered in |23] and shown to have a multifractal 
structure in the large N limit. For others values of 77, k the 
spectrum is in general not degenerate and this possibility does 
not exist. 

Notice that the eigenstates of —iF^ (with antiperi- 

odic boundarv conditions) in Eg . I34lare H\Q) andiJjl) where 
is the Hadamard operator 



H 



V2 U -1 



(36) 



Therefore, for this case A = H, Oq = Tr/2 and 9i = 0. The 
explicit expression of EQBM eigenstates with their respective 
eigenvalues for 3 qubits are: 



^1,0 
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with A: = 0,1,2. 



Notice that in this argument the role of the tensor product 



t (0.5,0.5) 



This 



of Hadamard gates is that of diagonalizing K 
is a peculiar case of the antisymmetric quantization (rj, k — 
1/2,1/2). Other values would still lead to a tensor product of 
one-qubit gates but with a more complicated structure given 
by Eq. |33 This is different from the use of the Hadamard 
transform in 1 6 ] where its role is to restore the parity symmetry 
to the shift map and is thus independent of the Floquet angles 
considered. 



B. Matrix of overlaps 

While the average PR gives a global measure of the similar- 
ity between the eigenbases of Bbvs (\^a)) and (l^'/s)) a 
more detailed view is given by the matrix of overlaps Af^a = 



8 



I (^'/sl^'a) |. If the eigenstates are ordered by their phases on 
the unit circle, the resulting matrix should be almost diago- 
nal with the off-diagonal elements signaling the importance 
of small components other than the main ones. 



EQBM's spectrum is not enough to accurately represent the 
s pectrum of QBMs. 




Figure 13: Matrix of overlaps (A/)^,^, = K^'/sl'&a)! with Floquet 
angles = k = 0.5 for i3 = 30 (left) and D = 32 (right). The 
eigenstates of both bases are ordered by growing eigenphases. 

Fig. ^] shows the matrix of overlaps for 13 = 30 and 
for the qubit case D — 32. For the latter the spectrum 
is concentrated at the values and therefore many 

eigenstates are degenerate. In this case we choose the basis 
as given in Eqs. |32]and|^ For I? = 30 on the other hand 
there are no degeneracies and the basis is unique. Fig. [T3l 
shows clearly that in both cases most states |<i>Q,) have a very 
large overlap with one basis state (this overlap being typically 
~ 0.9) while a few are mixtures of two. Some inversions in 
the ordering on the unit circle are also evident. 

The eigenfunctions are well approximated by the states that 
show no mixing in the matrix of overlaps. Fig. 1141 shows 
the comparison of three such eigenstates for 13 = 32 in the 
Husimi representation. The main difference in these eigen- 
functions are due to the lack of time reversal symmetry of 
the basis states. This symmetry, present in is reflected 
clearly in the pattern of zeroes which is symmetric with re- 
spect to the main diagonal. Notice that, while the basis states 
for D = 2^ can be analytically constructed, the basis of the 
EQBM for D ^ 2^ has to be found numerically. What we 
show here is that EQBM provides an excellent starting point 
to study the spectral properties of any baker quantization for 
all values of D. 



C. The spectrum 

Every QBM different from the EQBM has a complex spec- 
trum which follows the RMT predictions (non-degenerate, 
with eigenvalues repulsion, etc. |24|). In general, the spec- 
trum of EQBM will be "chaotic" too. But, when the Hilbert 
space is spanned by qubits and for special Floquet angles 
the EQBM has a short period and its spectrum is degenerate. 
Therefore, in these cases, the zeroth order approximation with 




Figure 14: Husimi representation for three eigenstates of the Bbvs 
(left) |<1>q) and *8 (l^'/j)) (right) in a Hilbert space spanned by 5 
qubits. 



The first order aproximation can be computed multiplying 
the zeroth order eigenvalues with the diagonal of the perturba- 
tion matrix (the diffraction kernel) represented in the EQBM 
basis. The corrected eigenvalues are then 



A 



(1) 



(38) 



We find that the diagonal elements of the diffraction kernel 
in Eq. |38]are close to the unit circle, their modulus tipically 
^ 0.9. We adopt their phase as first order correction to the 
A^^"* eigenphases. 

This approximation will be specially useful for the case of 
Hilbert spaces spanned by qubits since the corrections break 
degeneracy of the analytic EQBM's spectrum. Using Eas.l21l 
|22]and|^the first order approximation is 
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= ^ E e^'^'^'-KMS^'A^' 
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Figure 15: Eigenphases of the shift operator (S), the EQBM (EB), 
the first order correction of the EQBM with the Kernel B (EBpert ), 
and the BVS Baker (B) for dimension D = 32 (top) and D = 30 
(bottom). The continuous lines are the smoothed spectral densities. 



Kd.d^A'^^S"'W„) (39a) 

3=0 J 

I 

IV. CONCLUSION 

We have shown that many valid quantizations of the baker's 
map have a simple matrix (5B) as an essential core. This ma- 
trix has a structure that captures many features of the com- 
plexity of eigenfunctions of these maps. In spite of its simplic- 
ity, this matrix generally cannot be analytically diagonalized. 
The numerical eigenfunctions are excellent approximations to 
most baker map's eigenfunctions with an average participa- 
tion ratio of about 15 at Z? ~ 500. Many individual eigen- 
functions have participation ratios of values near 1. The eigen- 
values are also well approximated in their fluctuation proper- 
ties. The smoothed density of states is almost identical for all 
quantized families. The case £> = 2^ is very special. Some 
of its properties have been analyzed before as a toy model 1 13] 
or as the extreme case of the Schack and Caves family |20j. 
Its spectrum is highly degenerate and its semiclassical limit is 
different from all of the others. Both the eigenvalues and the 
eigenfunctions can be analytically obtained and the basis con- 
structed on primitive binary cycles. They have participation 
ratios that smoothly interpolate those of more general values 
of D. An important ingredient in S is the shift S. In fact, we 
could well have shifted the gate in Fig. |9]to the diffraction 
kernel and obtained a splitting in terms of the shift multiplied 
by a modified "perturbation". This procedure would yield as 
an "essential" baker just the shift S and would leave all the 
peculiarities of the different families to a modified diffraction 
kernel. The shift basis, parity projected, was recently adopted 
in ref. 01 • It has the advantage, over the EQBM basis, that it 
can be analytically constructed for all values of D using the 
number theory properties of the shift. However, as we show in 
Fig. [21 it consistently gives higher values participation ratio. 
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In Fig. El we compare the spectra obtained with this ap- 
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Appendix A: FACTORIZATION OF DFT WITH FLOQUET 
ANGLES 

In the text we repeatedly used a circuit representation 
for the Fourier transform. Here we provide the details of 
this representation. The discrete Fourier transform in a D- 



10 



dimensional Hilbert space with Floquet angles is 



this swap takes the permutation form 



1 



2ti 



{j\F'i^^\k) - ^ cxp[-z — (j + ry)(fc + n)] 



(Al) 



If D = D1D2, we can relabel the states as j = ji + Dij2 
and k 1)2/01 + with < ji,ki < Di - 1, and < 
j2, k2 < D2 — 1. Notice that index j2 is more significant than 
ji and viceversa for ki. We then obtain 



{j\F'D\Djk) - {.n\FZ\'\k^){j2\F^^^:\k2)eUi,k2) (A2) 

where, besides the Fourier transforms in each subspace, there 
is an an interaction term between ji and k2 given by 

The circuit of this factorization is represented in Fig. 







D1D2 










Figure 16: Circuit representation of the factorization of the quan- 
tum Fourier transform with Floquet angles in a Hilbert space with 
dimension D = D1D2 

The swap between the two subspaces is needed because of 
the different significance assigned to the Di, D2 subspaces in 
the decomposition of j and k. Asa matrix in the product space 



D-l 

SD1D2 = X! l-^i-^ 



D2 



{D-l)) 01 



(A4) 



where [x] is the integer part. Notice that Sd^Di — ^d\d2 
For Di = 2, S2D2 it is the shift operator utilized by Lakshmi- 
narayan |j8|]. 

Clearly the factorization is recursive so that it can be iter- 
ated until the different prime factors of D are reached. When 
D — 2^ and 77, k = the decomposition leads to the well 
known quantum circuit for the DFT I14ll in terms of Hadamard 
gates and diagonal phases. 

We give for reference the form of the interaction 8 for the 
special case of interest here of one qubit interacting with a 
subspace of dimension D {Di — 2, D2 — D) for the special 
values r],K — 0,0 and r/, n = 1/2, 1/2 

?7,K = 0,0 ^ e(e,fc) = e-'T^''' 
77, At =1/2, 1/2 ^ e(e,fc) =e*5e-*"('-^)(^-5) 

where e = 0, 1 and k = 0, . . . , D ~ 1. The first case gives the 
controlled phase familiar form the QFT, while the second is a 
more symmetrical interaction appropriate to the antisymmet- 
ric baker map. 

Using this decomposition the baker's map in the position 
representation can be split into the two kernels. This is done 
pictorially in Fig. ^\ The decomposition is used in Sec IID 
in the factorization of the QBM families. 
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Figure 17: Splitting of the BVS map in the product of two kernels. 
The thick line represents a space of dimension D/2 while the thin 
line is a qubit. The Floquet angles are implicit and the squares repre- 
sent the phase interaction in Eq IA3I 



